ADAPTIVE CONTROL SYSTEM HAVING DIRECT OUTPUT FEEDBACK 
AND RELATED APPARATUSES AND METHODS 

Cross-Reference to Related Application 

5 This application claims priority benefits of U. S. provisional application no. 

60/208,101 filed May 27, 2000 naming Anthony J. Calise, Naira Hovakimyan, and 
Hungu Lee as inventors. 

Statement of Government Rights in the Invention 
This invention was funded in part by the Air Force Office of Scientific 
10 Research (AFOSR) under Grant No. F4960-01- 1-0024. The United States 
Government therefore has certain rights in the invention. 

Field of the Invention 

Q The invention is directed to a system, apparatuses and methods for adaptively 

controlling a plant such as an aircraft, automobile, robot, or other controlled system. 
J P 15 Background of the Invention 

m Research in adaptive output feedback control of uncertain nonlinear dynamic 

JiJ systems is motivated by the many emerging applications that employ novel actuation 

: & devices for active control of flexible structures, fluid flows and combustion processes. 

Ifi These include such devices as piezoelectric films, and synthetic jets, which are 

m 20 typically nonlinearly coupled to the dynamics of the processes they are intended to 
:^ control. Modeling for these applications vary from having accurate low frequency 

models in the case of structural control problems, to having no reasonable set of 
model equations in the case of active control of flows and combustion processes. 
Regardless of the extent of the model accuracy that may be present, an important 
25 aspect in any control design is the effect of parametric uncertainty and unmodeled 
dynamics. While it can be said the issue of parametric uncertainty is addressed within 
the context of adaptive control, very little can be said regarding robustness of the 
adaptive process to unmodeled internal process dynamics. 

Synthesis approaches to adaptive output feedback control typically make use 
30 of state estimation, and therefore require that the dimension of the plant is known. 
Some approaches further restrict the output to have full relative degree, or restrict the 
uncertainties in the plant to be an unknown function of the output variables. It would 
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be desirable to remove all these restrictions by adopting a direct output feedback 
approach that does not rely on state estimation. One of the immediate consequences 
of such an approach would be that the dimension of the controlled plant need not be 
5 known. Consequently, the resulting system would be applicable to plants having both 
parametric uncertainty and unmodeled dynamics. Furthermore, it would be desirable 
to produce a control system that is not only robust to unmodeled dynamics, but also 
learns to interact with and control these dynamics. 

Output feedback control of full relative degree systems was introduced by 

10 Esfandiari and Khalil, 1992, "Output feedback stabilization of fully linearizable 
systems," InternationalJournal of Control, 56(5): 1007-1037. In their publication the 
authors formulated a control methodology that involves a high gain observer for the 
reconstruction of the unavailable states. A solution to the output feedback stabilization 
problem for systems in which nonlinearities depend only upon the available 

15 measurement, was given by Praly, L. and Jiang, Z. (1993)/ "Stabilization by output 
feedback for systems with iss inverse dynamics," System & Control Letters , 21:19- 
33. Krstic, M., Kanellakopoulos, I., and Kokotovic, P. (1995), Nonlinear and 
Adaptive Control Design, John Wiley & Sons, Inc. New York and Marino, R. and 
Tomei, P. (1995). Nonlinear Control Design: Geometric, Adaptive, & Robust, 

20 Prentice Hall, Inc., presented backstepping-based approaches to adaptive output 
feedback control of uncertain systems, linear with respect to unknown parameters. An 
extension of these methods due to Jiang can be found in Jiang, Z. (1999), A combined 
backstepping and small-gain approach to adaptive output feedback control. 
Automatical 35:1131-1139. 

25 For adaptive observer design, the condition of linear dependence upon 

unknown parameters has been relaxed by introducing a neural network (NN) in the 
observer structure of Kim, Y. and Lewis, F. (1998), High Level Feedback Control 
with Neural Networks, World Scientific, NJ. Adaptive output feedback control using a 
high gain observer and radial basis function neural networks (NNs) has also been 

30 proposed by Seshagiri, S. and Khalil, H. (2000), "Output feedback control of 
nonlinear systems using {RBF} neural networks," IEEE Transactions on Neural 
Networks, ll(l):69-79 for nonlinear systems, represented by input-output models. 
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Another method that involves design of an adaptive observer using Amotion 
approximators and backstepping control can be found in Choi, J. and Farrell, J. 
(2000), "Observer-based backstepping control using on-line approximation," 
Proceedings of the American Control Conference, pages 3646-3650. However, this 
result is limited to systems that can be transformed to output feedback form, i.e., in 
which nonlinearities depend upon measurement only. 

The state estimation based adaptive output feedback control design procedure 
in the Kim and Lewis 1998 publication is developed for systems of the form: 

x =f(x) + g(x)£ c (1) 
y = x dim x = dim y = dim u, (2) 

which implies that the relative degree of y is 2. In Hovakimyan, N., Nardi, F., Calise, 
15 A., and Lee, H. (1999), "Adaptive output feedback control of £Lclass of nonlinear 
systems using neural networks," International Journal of Control that methodology is 
extended to full vector relative degree MIMO systems, non-affine in control, 
assuming each of the outputs has relative degree less or equal to 2: 

20 x=f(x 9 S c ) (3) 

y = h(x) dim y = dim u < dim x. (4) 

These restrictions are related to the form of the observer used in the design procedure. 
25 Constructing a suitable observer for a highly nonlinear and uncertain plant is not an 
obvious task in general. Therefore, a solution to adaptive output feedback control 
problem that avoids state estimation is highly desirable. 

Summary of the Invention 
The adaptive control system (ACS) and method of this invention uses direct 
30 adaptive output feedback to control a plant. The system can comprise a linear 
controller (LC) and an adaptive element (AE). The linear controller can be used as a 
dynamic compensator to stabilize a model of the plant, and provide output regulation. 
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The adaptive element can compensate for disturbances, and modeling error resulting 
from approximation in modeling of the plant. The adaptive element can comprise a 
neural network (NN). The adaptive element can receive a signal from the linear 
5 controller used to adapt its NN's weights. The input vector to the NN. can comprise 
current and/or past plant output signals together with other available signals. The past 
plant output signal(s) can be used as inputs to the NN to ensure boundedness of the 
adaptive element in controlling the plant. The adaptive control system can comprise 
an error conditioning element having a low-pass filter designed to satisfy a strictly 

10 positive real (SPR) condition of a transfer function associated with Lyapunov stability 
analysis of the control system. The stability analysis can be used to construct the NN 
adaptation law using only the plant output signal(s) and other available signals as 
inputs to the NN, and to ensure boundedness of error signal(s) of the closed-loop 
adaptive control system. Apparatuses forming components of the ACS are also 

15 disclosed. ' ' 

A method of the invention comprises generating at least one control signal 8 C 
to regulate a plant output signal y by feedback of the plant output signal y, and 
optionally other sensed variables related to the state of the plant, in which y is a 
function of the plant state having known but unrestricted relative degree r. The control 

20 signal S c can be generated so as to control the plant based on an approximate 
dynamic model, and so as to control the plant in the presence of unmodeled dynamics 
in the plant based on an adaptive control technique. The adaptive control technique 
can be implemented with a neural network. Related methods are also disclosed. 

These together with other objects and advantages, which will become 

25 subsequently apparent, reside in the details of construction and operation of the 
invented system, methods, and apparatuses as more fully hereinafter described and 
claimed, reference being made to the accompanying drawings, forming a part hereof, 
wherein like numerals refer to like parts throughout the several views. 

Brief Description of the Drawings 

30 Fig. 1 is a general block diagram of an adaptive control system for controlling 

a plant based on a plant output signal y in accordance with the invention; 
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Fig. 2 is a schematic diagram of the adaptive control system reduced to 
elements relevant to tracking error signal analysis; 

Fig. 3 is a relatively detailed view of a linear dynamic compensator of the 
5 adaptive control system; . . . _ 

Fig. 4 is a relatively detailed view of an adaptive element; 
Fig. 5 is a block diagram of a processor-based adaptive control system using 
direct output feedback; 

Fig. 6 is a flowchart of a general method of the invention; 
10 Fig. 7 is a graph of commanded output signal y c (broken line) and plant output 

signal y (solid line) versus time for a control system without unmodelled dynamics 
using only a linear compensator to control a plant; 

Fig. 8A is a graph of commanded output signal y c (broken line) and plant 
output signal y (solid line) versus time for a control system without unmodelled 
15 dynamics and with the adaptive element; - - ^ 

Fig. 8B is a graph of adaptive control signal v ad and inversion error signal 
A for a control system without unmodelled dynamics and with the adaptive element; 

Fig. 9 is a graph of commanded output signal y c (broken line) and plant output 
signal y (solid line) versus time for a control system with unmodelled dynamics and 
20 with the adaptive element; 

Fig. 1 OA is a graph of commanded output signal y c (broken line) and plant 
output signal y (solid line) versus time for a control system with unmodelled 
dynamics, and with the adaptive element; and with the linear controller. 

Fig. 1 OB is a graph of commanded output signal y c (heavy line), plant output 
25 signal y (line with relatively moderate oscillations), and connection weights W (line 
with relatively heavy oscillations) versus time for a control system with unmodelled 
dynamics and with an adaptive neural network element and linear controller. 

Detailed Description of the Invention 
As used herein, the following terms have the following definitions: 
30 "Actuator" can be virtually any device capable of affecting the state of a plant 

to control one or more degrees of freedom thereof. Such actuator can be a motor, 
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motor-driven screw, a hydraulic cylinder, a pump or valve controlling a stream of air, 
a thermal heater, a compressor or suction generator, or other device. 

"Adaptive control system" means a control system having the capability to 
5 adapt to changes in a controlled plant or its environment over time. 

"And/or" means either one or both of the elements before and after this term. 
Hence, "A and/or B M means "A" or "B" or "A and B". 

"Direct output feedback" refers to a control system, apparatus or method that 
employs feedback of an "output" that is a function of the full state(s) existing in a 
10 plant. "Direct" refers to the fact that no state estimation is used to estimate plant states 
not present in or not derivable from the "output". 

"Operator" can be a human or computer, that receives and input and generates 
and output based on the current and past history of the input., for example, senses a 
plant output using a plant output signal, and generates a commanded state signal to 
15 control the plant. * ^ 

"Memory" can be a random-access memory (RAM), read-only memory 
(ROM), erasable read-only programmable memory (EPROM), or other memory 
device capable of storing a control program and data executable by a processor. 

"Plant" refers to a system controlled by a control system. For example, the 
20 plant can be an aircraft, spacecraft, space-launch vehicle, satellite, missile, guided 
munition, automobile, or other vehicle. The plant can also be a robot, or a pointing or 
orientation system such as a satellite orientation system to orient power-generation 
panels, a transceiver, or a docking mechanism. Such plant can also be a braking 
system, an engine, a transmission, or an active suspension, or other vehicle subsystem. 
25 The plant can be a manufacturing facility or a power generation facility. In general, 
the plant could be virtually any controllable system. 

"Processor" can be a microprocessor such as a Xeon® or Pentium® brand 
microprocessor produced by Intel® Corporation, an Athlon® brand microprocessor 
commercially available from AMD® Corporation, Sunnyvale, California, which can 
30 operate at one (1) megahertz or more, a microcontroller, a field programmable gate 
array ("FPGA"), a programmable logic array ("PLA"), a programmed array logic 
("PAL"), or other type of data processing or computing device. 
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"Relative degree" applies to a regulated variable (such as plant output signal y) 
and corresponds to the number of times the variable must be differentiated with 
respect to time before an explicit dependence on the control variable (such as the 
5 command control signal 5 C ) is revealed. 

"Sensor" can be virtually any device(s) for sensing a degree of freedom of a 
plant's state, whether alone or in combination with one or more other sensors. The 
sensor can be virtually any device suitable for sensing information regarding a plant's 
state. For example, the sensor could be a gyroscope for detecting orientation of a 
10 vehicle such as an aircraft, i.e., pitch or roll attitudes or side slip. The sensor can also 
be a temperature or pressure sensor, a position, velocity, or inertial sensor. 

"(s)" means one or more of the thing meant by the word preceding "(s)". Thus, 
basis function(s) means one or more basis functions. 

"State" refers to a property of a plant to be controlled which is sufficient to 
15 completely define the condition of the plant at any time instant. For example, 
elements of the state can be a position, velocity, acceleration, mass, energy, 
temperature, pressure, volume, etc. of an object associated with a plant that is to be 
controlled. 

"State feedback" pertains to a situation in which the entire state of the plant 
20 can be sensed and used to control the plant through feedback. 

"Strictly positive real" is a property that pertains to the transfer function of a 
linear time-invariant system. The transfer function, G(s), is a ratio of polynomials in 
the variable, 's', which is a complex variable having a real and imaginary part. Let s = 
a + jb, were a is the real part and b is the imaginary part. Then the transfer function is 
25 called 'Strictly Positive Real 1 if the following two conditions are satisfied: 

1) G(s) is asymptotically stable (all the poles of G(s) have real parts < 0); and 

2) G(jb) > 0 for all values of the real variable V. This definition can be found in 
Khalil, H.K., "Nonlinear Systems, Second Edition, Prentice-Hall, 1996, p. 404. 

"Variable" refers to any signal that can be changed independently of the plant 
30 states, such as the control variable, or that dependent upon time either directly, or 
indirectly because it depends upon plant states that are time varying, such as the 
output variable. 




1. General Description of Adaptive Control System and Method 

As shown in Fig. 1, an adaptive control system (ACS) 10 can be used to 
5 control ..a plant 12 using 'direct output feedback 1 as opposed to 'state feedback/ In 
contrast, the disclosed adaptive control system 10 uses 'direct output feedback 1 in 
which a function of the full plant state, as opposed to all plant states or estimates 
thereof, is used for feedback control of the plant 12. The adaptive control system 10 of 
this invention is therefore highly useful in the control of plants, especially non-linear 

10 plants in which the full state of the plant cannot be sensed through practical or 
economically-feasible devices or techniques, particularly if the plant contains 
unmodeled or unknown dynamics. 

In Fig. 1, the ACS 10 comprises a linear controller (LC) 14, stable adaptive 
element (AE) 16, and model inversion unit (MIU) 18. The ACS 10 can further 

15 comprise error signal generator (ESG) 20 and summing unit 22^ The ACS 10 can 
further comprise an operator 20, operator interface unit 22, and command filter unit 
24. These elements permit the operator 22, which can be human or machine, to 
interact with ACS 10 to control the plant 18. The plant 12 comprises a system 30 to be 
controlled. In general, the system 30 is a non-linear system, although it can be linear. 

20 The linear controller 14 is implemented to affect approximate control of the plant, 
optionally using linear control. The AE 16 serves to implement adaptive control of 
nonlinearities of the plant 12 that are not compensated by the linear controller 14. The 
combined effect of the LC 14 and AE 16 is used to control the plant 12. 

The ACS 10 is now described in further detail. The operator interface unit 26 

25 is coupled to receive the plant output signal y which, as previously described, is a 
function of the full state of plant 12 having known but unrestricted relative degree r in 
which r is the number of times the plant output signal y must be differentiated in order 
for the corresponding control variable, 8 C , to become explicit. The operator interface 
unit 26 can be an electronic interface between an ACS bus and the operator 24 if 

30 implemented as a processor, for example, or can be a display, gauge, meter, light, or 
other indicator if the operator 24 is human. The operator 24 generates command 
signals based on the plant output signal y from the operator interface unit 26. The 
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operator 24 generates a command action or command signal based on the plant output 
signal y. The operator 24 supplies the command action or command signal to the 
command filter unit 28 that generates filtered commanded signals y c and y { c r) based 
5 thereon. The integer r denotes the relative degree of the regulated variable, and signal 
y { c r) denotes the derivative with respect to time of the filtered commanded signal 
y c . The command filter unit 28 is coupled to supply the filtered commanded signal 
y c to the error signal generator 20. The error signal generator 20 is also coupled to 
receive the plant output signal y. Based on the filtered commanded signal y c and the 
10 plant output signal^, the error signal generator 20 generates a tracking error signal y . 
The error signal generator 20 is coupled to supply the tracking error signal y to the 
linear controller 14. 

The linear controller 14 generates a pseudo-control component signal v dc 
based on the tracking error signal y by operating on such error signal with a summing 

15 node, feedback network, and gain multipliers implementing the transfer function 
Ndc( s ) / D dc( s )- The linear controller 14 also generates a transformed signal y ad 
based on the tracking error signal y by operating on such error signal with a summing 
node, one or more integrators, and a feedback network from output terminal(s) of the 
integrators) that have gain multipliers) implementing the transfer function N a d(s) / 

20 Dd c (s). The linear controller 14 generates the transformed signal y ad so that the 
transfer function from the transformed signal y ad to the adaptive control signal v a d is 
strictly positive real (SPR). The linear controller 14 is coupled to supply the 
transformed signal y ad to the AE 16. 

More specifically, the error conditioning element 38 of the AE 16 is coupled to 

25 receive the transformed signal y ad . The error conditioning element 38 is also coupled 
to receive basis function(s) <J> and generates the training signal 5 based on the signal 
y ad and the basis function(s). The error conditioning element 38 can generate the 
training signal 8 by filtering the basis function(s) <|) and multiplying the resulting 
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signal <|> f by the transformed signal y ad . The error conditioning element 38js coupled 
to supply the training signal 5 to the neural network adaptive element (NNAE) 36. 

The NNAE 36 uses the training signal 5 to adjust connection weights W of its 
5 neural network to adapt to plant dynamics that are unmodeled and therefore not 
adapted to by the LC 14. The delay element 40 of the AE 16 is coupled to receive the 
pseudo-control signal v and the plant output signal y and is coupled to supply non- 
delayed signals v, y, and delayed versions v^, yd thereof generated by the delay 
element 40 based on the signals v, y to the NNAE 36. The delayed signals yd are 

10 delayed relative to the plant output signal^ by time delay increments d to {n - 1) • d 9 n 
being the number of the full plant states, although in practice fewer or more such 
delays may be used. The delayed signal(s) vfi are delayed relative to the pseudo- 
control signal v by time delay increments d to (n - r - I) - d , r being the relative 
degree. The use of these delayed signals assures that the weight coefficient(s), W, and 

15 error signal(s) y are bounded so that the ACS 10 is stable. 

Based on the connection weight(s) W determined by training signal 6, and the 
pseudo-control signal v and delayed version(s) vj thereof and/or the plant output 
signal y and delayed version(s) yd thereof supplied as inputs to the NNAE 36, the 
NNAE generates the adaptive control signal v a d- The NNAE 36 is coupled to supply 

20 the adaptive control signal v a d to the summing unit 22. The summing unit 22 is also 
coupled to receive the pseudo-control component signal vdc fr° m ^ linear controller 
14 and the r-th time derivative of the commanded state signal y { c r) from the command 
filter unit 28. Based on the signals y ( c r) , vdc > v ad > the summing unit 22 generates the 
pseudo-control signal v. The summing unit 22 is coupled to supply the pseudo-control 

25 signal v to the model inversion unit 18. 

The model inversion unit 18 is also coupled to receive the plant output signal 
y. The model inversion unit 18 generates a command control signal 8 C based on the 
pseudo-control signal v and the plant output signal y. More specifically, the model 
inversion unit 18 subjects the signals v, y to a function that inverts the linear control 

30 model implemented by the linear controller 14, to produce the command control 
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signal 8 C . The model inversion unit 1 8 is coupled to supply the resulting^ command 
control signal 8 c to the actuator(s) 32 of the plant 12. The actuator(s) 32 are coupled to 
or associated with the controlled nonlinear system 30 so that the control the state(s) of 
5 such system, based on the command control signal 8 C . The sensor(s) 34 are coupled or 
associated with the controlled nonlinear system 30, and generates the plant output 
signal y that is a function of the full plant state(s) of the controlled nonlinear system 
30. 

Model inversion in the unit 18 can be performed in the following way. For the 
10 scalar case (p=l) if the regulated output, y 9 has relative degree r , then the output 
equation can be differentiated r times with respect to time until the control appears 
explicitly. In this case, we assume that we have a model of this derivative 

y (r) =h r (y,5 c ) = v . (5) 

15 

Thus the model inversion of the unit 18 is defined by 

s c =K\y>v). (6) 

20 2. Specific Description of Adaptive Control System and Method 

Let the dynamics of an observable nonlinear single-input-single-output (SISO) 
system be given by the following equations: 

x=f(x,J c ), y-h(x) (7) 

25 

where x e Q c 9t n is the state of the system, 5 c ,y e 91 are the system input (control) 
and output (measurement) signals, respectively, and f(y),h(')eC w are unknown 
functions. Moreover, n need not be known. 

To ensure proper operation, the following assumption is made in the 
3 0 implementation o f the ACS 1 0 . 



11 



Assumption 1. The dynamical system of Eq. (7) satisfies the output feedback 
linearization conditions Isidori, A. (1995), Nonlinear Control Systems, Springer- 
Verlag, Inc. with relative degree r, i.e., 

5 

y (r) =h r (x,S c ). (8) 

Here, h A—, such that ^ = 0 for 0 < i < r and^- * 0. 
r== dt r du du 

10 This disclosure addresses the design of an output feedback control law that utilizes the 
available measurement y(t), to obtain system output tracking of a bounded trajectory 
y c (t) that is assumed to be r-times differentiable, i.e., y c e C r . The difference between 

unknown dynamics function 1^ and its estimate h r , or the modeling error, is mapped 
using a NN. This mapping has to be based on measured input and output data only. 

15 To this end, the universal approximation property of neural networks and the 
observability of the system are utilized to construct this mapping on-line using 
measured input/output time histories. These various features of the proposed control 
design scheme are presented in the next section. 

3 Controller Design 

20 3.1 Feedback Linearization 

Feedback linearization is performed by introducing the transformation 

v = h r (y,8 c ), (9) 

25 where v is commonly referred to as a pseudo control signal, and h r (y,8 c ) is the best 
available approximation of h r (x,8 c ). Then, the system dynamics can be expressed 
as 

y (r) =v + A', (10) 

30 where 
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A'-A'(x,5 c ) = h r (x,5 c )-h r (y,5 c ). m (11) 

Using this transformation, Eq. (10) represents the dynamic relation of r integrators 
5 between the pseudo-control signal v and the plant output signal y, where the error A' 
acts as a disturbance signal. Assuming that the plant output signal y is required to 
track a known bounded input command signal y c9 the pseudo-control signal v is 
chosen to have the form 

10 v = y< r) +v dc -v ad , (12) 

where v dc is the output of a stabilizing linear dynamic compensator for the linearized 
dynamics in Eq. (10) with A' = 0, and is the adaptive control signal designed to 
approximately cancel A' . The r-th derivative of the input signal, y< r) , is introduced as 
15 a feedforward term to condition the error dynamics. This derivative can be easily 
obtained if the tracking (or command) signal^ is generated using an >-th (or higher) 
order reference model forced by an external input. The reference model serves to 
define the desired response of the closed loop system. The input to the dynamic 
compensator is the tracking error, which is defined by 
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y = y c -y. 03) 



It is important to point out that the model approximation function h r ( v ) should be 
defined so that it is invertible with respect to u, allowing the actual control input to be 
25 computed by 

5 c =hr 1 (y,v). (14) 
Clearly, the accuracy of the approximation h r (x,h^(y,v))- v is governed by 

30 
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A'(x,8 c ) = A(x,y,v) = h r (x,h r ^y, v))- v . 



(15) 



From Eq. (11) and Eq. (15), notice that A depends on v ad through v, whereas v ad has to 
be designed to cancel A. The following assumption is introduced to guarantee 
existence and uniqueness of a solution for v^: 

10 Assumption 2. The map h-> A is a contraction over the entire input domain of 
interest. 

Using Eq. (1 1), the condition in Assumption 2 implies: 
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dA 



dv 



ad 



5(h r -h r ) 55 c dv 



du dv dv a d 



5(h r -h r )5S c 



du dh T 



<1 



(16) 



which can be re-written in the following way: 



dh r /db c 



dh T /d8 c 



-1 



<1 
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The condition (17) is equivalent to the following two conditions 



(17) 



(1) sgn(dh r /ao c )=sgn(dh r /6S c ) 

(2) Idh r /d8 c >|dh r /S8 c |/2>0. 



25 The first condition means that control reversal is not permitted, and the second 
condition places a lower bound on the estimate of the control effectiveness in (14). 
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3.2 Control System Architecture 

Based on the above description, the overall control system architecture is 
5 presented in Fig. 1. The central components of the system are: (a) the model 

inversion/linearization unit 18 implementing block h~ l (y 9 8 c ) 9 (b) the adaptive neural 
network based element 16 is designed to minimize the effect of A, and (c) the linear 
dynamic compensator of the linear controller 14. The input into the ACS 10 is the 
reference command tracking signal y c and its r-ih derivative y< r) , generated by, e.g., a 

10 reference model forced by an external input. 

It is important to note the two output signals (v dc ,y ad ) of the linear 
compensator. The pseudo-control component signal, v dc , is designed to stabilize the 
linearized system, as described earlier. The transformed signal, y ad , is a linear 
combination of the compensator states and its input, i.e., the tracking error signal y . 

15 This signal is generated to ensure an implementable error signal 8 that is used to adapt 
the NN weights W of the NNAE 36. 

3.3 Tracking Error Signal Analysis 
The analysis presented in this subsection is carried out to facilitate the design 
of the NNAE 36 and the second output signal y a d of the l inear dynamic compensator 

20 14. To formulate the overall tracking error dynamics of the controlled system, the 
specific choice of the pseudo-control signal v is given by Eq. (12) is substituted into 
Eq. (16), leading to 

y^y^+v^-v^ + A, (18) 

25 or alternatively 

y (r) =-v, c +v a ,-A. (19) 

These error dynamics are depicted schematically in Fig 2. More specifically, under 
tracking error dynamics analysis, the ACS 10 reduces to a summing node 42, an 
30 integrator 44, and a linear dynamic compensator (LDC) 46. The summing node 42 is 
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coupled to receive the signal v ad - A and the pseudo-control component signal v^. The 
summing node 42 subtracts the signal from the signal v ad - A to generate the signal 
y (r) . The r-th degree integrator 44 integrates the signal y (r) to produce the signal y . 
5 The LDC 46 is coupled to receive the signal y from the integrator 44. Based on the 
signal y , the LDC 46 generates the signal v dc that is fedback to the summing node 42. 
The LDC 46 also generates the signal y ad based on the signal y. The LDC 46 
implements a transfer function that is SPR to map the signal y to the signal v ad - A to 
ensure stability of the ADC 10. 



where s represents the complex Laplace variable. The LDC 46 can comprise transfer 
function elements 48, 50. The transfer function element 48 can be used to implement 
the transfer function N dc (s)/D dc (s) mapping the signal y to the signal v dc . The 
transfer function element 50 can be used to implement the transfer function 
20 N ad (s) I D dc (s) map the signal y to the signal y ad . Further details regarding the LDC 
46 are described below. 

Assumption 3. The linearized system in Fig. 2 is stabilized using a stable linear 
dynamic compensator 46, i.e., the roots of the denominator polynomial Ddc( s ) are 
25 located in the open left half plane of the complex plane s. 

Since the linearized system dynamics, and hence the error dynamics, consist of r pure 
integrators, this assumption introduces only a very mild restriction on the design. 
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The single-input two-output transfer matrix of the linear dynamic compensator is 
denoted by 




(20) 



15 



16 



Based on the compensator defined in Eq. (20), the closed loop transfer function of the 
system depicted in Fig. 2 is given by: 



y a d(s)=— 



-(v ad - A)(s) AG(s)(v ad - A)(s) . 



s'D dc (s) + N dc (s) 



(21) 



Analyzing the denominator of Eq. (21), the Routh-Hurwitz stability criterion implies 
that a necessary condition for closed loop system stability is that the degree of the 
compensator numerator, N dc (s), and hence of its denominator, D dc (s) , should be at 
10 least (r - 1), i.e., 



qAdeg(D dc (s)) > deg(N dc (s)) > r - 1 . 



(22) 
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This dictates the design of the linear dynamic compensator 



v dc (s) = ^^y(s), 



(23) 



which can be carried out using any linear control design technique (classical, pole 
placement, optimal LQ, etc.), with the constraint of assumption 3. The numerator 
20 N ad (s) , associated with the output y ad , does not affect the stability of the error system 
of Fig. 2. 

3.4 Neural Network Based Approximation 

The term "artificial neural network" has come to mean any architecture that 
25 has massively parallel interconnections of simple "neural" processors. Given x e R Nl , 
a three layer-layer NN has an output given by: 



N 2 



, i = l,...,N, 



(24) 
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where <(>(.) is the activation function, v jk are the first-to-second layer interconnection 
weights, and w- are the second-to-third layer interconnection weights. 9^ and 0^ 
are bias terms. Such an architecture is known to be a universal approximator of 
5 continuous nonlinearities with squashing activation functions. See Funahashi, K. 
(1989), On the approximate realization of continuous mappings by neural networks. 
Neural Networks, 2:183-192; Hornik, K., Stinchcombe, M, and White, H. (1989), 
Multilayer feedforward networks are universal approximators, Neural Networks, 
2:359-366. 

1 0 Linearly parameterized neural networks 

y = W T 0(x) (25) 

are universal approximators as well, if vector function 0(*) can be selected as a basis 
15 over the domain of approximation. Then a general function f (x) e C k , 
xeDc5R n can be written as 

f(x) = W T 0(x) + s(x), (26) 

20 where s(x) is the functional reconstruction error. In general, given a constant real 
number 8* > 0, f(x) is within s* range of the NN, if there exist constant weights W, 
such that for all xeiR n Eq.(20) holds with ||s|| < e* . 

Definition 1. The functional range of NNAE 36 is dense over a compact domain 
x € D , if for any f (•) e C k and e* there exists a finite set of bounded weights W, such 

25 that Eq. (26) holds with ||e|| < e* . 

Various publications show that the functional range of NN in Eq. (25) is dense 
for different activation functions 0(-). See Cybenko, G. (1989) publication. 
Approximation by superpositions of sigmoidal function, Mathematics of Control, 
Signals, Systems, 2(4):303-314; Park, J. and Sandberg, I. (1991), Universal 

30 approximation using radial basis function networks, Neural Computation, 3:246-257 r ; 
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Sanner, R. and Slotine, J. (1992), Gaussian networks for direct adaptive control, IEEE 
Transactions on Neural Networks, 3(6):837-864. 

The following theorem extends these results to map the unknown dynamics of 
5 an observable plant from available input/output history. 

Theorem 1. Given e* >0, there exists a set of bounded weights W, such that 
A(x,y,v), associated with the system (l)-(5), can be approximated over a compact 
domain DcQxR by a linearly parameterized neural network 

10 A = W T 0(r|) + eO 1 ), ||s||<e* (27) 

using the input vector 

ri(t) = [l v d T (t) yl(t)f, (28) 

15 where 

vj (t) = [v(t) v(t - d) • • • v(t - (n, - r - l)d) J 

yj(t) = [y(t) y(t-d) - y (t-( ni -i)d)f 

20 with n, > n and d > 0 , provided there exists a suitable basis of activation functions 
#(•) on the compact domain D. 

The output of the adaptive element 16 in Fig. 1 is designed as 

v ad =W T 0(r|), (29) 

25 

where W are the estimates of the weights. Eq. (29) will always have at least one fixed- 
point solution, so long as 0(-) is made up of bounded basis functions. 



30 
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3.5 Construction of SPR Transfer Function 



5 



As discussed earlier, the second output of the linear dynamic compensator 46, 
y ad , will be used to construct the rule for adapting W in Eq. (29). Using Eqs. (27) and 
(29) in Eq. (21) implies: 



y ad (s) = G(s)(w T 0(ti)-6) 



(30) 



where W = W - W is the weight error. As will be seen in the next section, for the NN 
10 adaptation rule to be realizable, i.e. dependent on available data only, the transfer 
function G(s) must be strictly positive real (SPR). However, the relative degree of 
G(s) is at least r. When the relative degree of G(s) is one, it can be made SPR by a 
proper construction of N ad (s). If r >1, G(s) cannot be SPR through this technique 
alone. 

15 To achieve SPR in the r >1 case, following the Kim and Lewis, 1998 

publication, a stable low pass filter T _1 (s) is introduced in Eq. (30) as:' 



20 where 0 f and e f are the signals 0 and 8, respectively, after being filtered 
through T" 1 (s) , and 8 m (s) is the "so-called" mismatch term given by 




(31) 




(32) 
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that can be bounded as 



||8 m (t)||<c|w| F) c>0. 



(33) 
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The numerator of the transfer function G(s)T(s) = G(s) in Eq. (31) is J(s)N ad (s) . 
The polynomial T(s) is Hurwitz, but otherwise can be freely chosen, along with the 
numerator polynomial N ad (s)of Eq. (21) that defines the output y ad . Hence, the 
5 numerator polynomial (or the zeros) of G(s) can be freely chosen to make it SPR. 

Two approaches can be utilized in constructing T(s) and N ad (s) to makeG(s) SPR. 
To avoid an unrealizable feedthrough, G(s) will be assigned (r + q - 1) zeros, thus 
making it relative degree one. 



10 Zero placement approach: Since G(s) is a stable transfer function, its zeros can be 
easily placed to make it SPR, e.g., by interlacing them with its poles. From Bode plot 
analysis it is easy to conclude that such a pole-zero pattern will ensure a phase shift in 

the range of ± 90° . 



15 LKY approach: Assume that 



_ b^ -1 +b 2 s p - 2 +--- + b 

G(s) = -4 |u L (34) 

s p +a,s p +-» + a p 



where p = r + q is the number of the closed loop poles. The controller canonical state 
space realization of this transfer function is given by 

20 z = A cl z + B cl (w T 0 f + 6 - s f ) 

y ad =C cl z, (35) 

where 





"-a, -a 2 ... -a p 




V 




1 0 — 0 




0 


A cl = 


0 '•. "-. : 


B d = 


0 




0 ••• 0 1 0 




0 
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C cl =[ b, b 2 ••• b p J 

G(s) is SPR if and only if it complies with the Lefschetz-Kalman-Yakubovitz (LKY) 
5 Lemma, Ioannou, P. A. and Sim, J. (1996), Robust Adaptive Control, Prentice Hall, 
Inc., p. 129, i.e., there exists Q>-0 such that the solution P of 

A2iP + PA d =-Q (36) 

is positive definite and 

10 

PB cl =Cj,. (37) 

G(s) can be constructed utilizing the LKY condition as follows: 

1 5 a) Choose Q 0 and solve Eq. (36) for P >- 0 . 

b) Using Eq. (37), compute C cl , which in this canonical form is 

simply the first column of P. From Eqs. (33) and (34), the elements of 
C cl are also the coefficients of the numerator polynomial ofG(s) . Since 

G(s) is SPR, it is guaranteed that this numerator is Hurwitz. 
20 c) Solve the numerator polynomial for its roots. 

From the zeros obtained by either of the above methods, choose (r - 1) of these to 
construct T(s), while the remaining q zeros makeup N ad (s). The fact that the 
numerator of G(s) is Hurwitz ensures also that T(s) and N ad (s) are individually 
25 Hurwitz. There is freedom in scaling T(s) and N ad (s), which could be utilized to 
normalize the maximum gain of T" 1 (s) . 

To summarize, N dc (s)/D dc (s) is designed to stabilize the linearized system 
dynamics, while N ad (s)is constructed to meet the SPR condition needed for a 
realizable implementation. 
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Neural Network Adaptation Rule 

5 As is evident from Eq. (3 1 ), the filter T" 1 (s) should operate on all the components of 
the NN vector 0 . All these filters can be cast in one state space realization: 

. z f = A f z f +B f 0 

(38) 

10 0 f =C f z f , 

where the diagonal blocks of the state space matrices (A f ,B f ,C f ) are constructed 
from a state space realization of the filter T _1 (s) . Since the filter is stable, 3 P f >- 0 , 
satisfying , . 

15 

AjP f +P f A f =-Q f (39) 

for any positive definite Q f >- 0 

20 The signals 0 f are used in the following NN W weight adaptation rule 

dW/dt = -F|y ad 0 f + A, w w], (40) 

where F >- 0 and X w > 0 are the adaptation gains. In the next section it is proven that 
25 this adaptation rule ensures boundedness of the system error signals and the network 
weights. The NNAE 16 of Fig. 1 is depicted in more detail in Fig. 4. 

Fig. 4 is an exemplary embodiment of the adaptive element 16 provided by 
way of example and not limitation as to possible implementations of the NNAE 16. 
The time delay element 40 comprises one or more time-delay (TDL) units 52 1 -52 nl _ r _ 1 
30 coupled to receive the pseudo-control signal v, and TDL units 54,-54 nM coupled to 
receive the plant output signal y. The TDL units 52 r 52 nl _ r _i, 54,-54 nM generate delayed 
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versions v d , y d of the signals v, y , and are coupled to supply these delayed signals 
v d , y d as well as undelayed signals v, y 9 to the NNAE 36. The neural network (NN) 
64 of the NNAE 36 multiplies the signals v, v dJ y 9 y d by respective weight data V and 
5 transmits the resulting signals to respective basis functions <)>(•) 56 l5 56 2 ,..., 56 NI . 
The basis functions (j>(-) 56, , 56 2 ,..., 56 NJ are coupled to receive V-weighted signals 
v, v d , y 9 y d and generate respective signals based thereon. The generated signals are 
multiplied by respective weight data W and summed at respective summation nodes 
58„ 58 2 , 58 N2 . The NNAE 36 is coupled to supply the resulting summed signals as 
10 the vector signal v ad to the summing node 22 of Fig. 1 for generation of the pseudo- 
control signal v. 

To ensure boundedness of the basis functions <)>(•) 56,, 56 2 ,.. 56 N1 and neural 
network weights V,W, the NNAE 36 is coupled to supply the basis functions ((>(•) 56„ 
56 2 ,..., 56 N1 as signals to the error conditioning element 3S : The error conditioning 

15 element 38 comprises a filter 60 and a multiplier 62. The filter 60 operates on the 
basis functions 56 1? 56 2 ,...,56 N1 with a filtering transfer function T l (s) as previously 
described with respect to Equation (32) to generate filtered basis functions ty { (•). The 
filter 60 is coupled to supply the filtered basis functions <|> f (•) to the multiplier 62. The 
multiplier 62 is also coupled to receive the transformed signal j> ad . The multiplier 62 

20 generates the signal 5 that is a vector product of the signals <|> f (•)> J> ad . The multiplier 62 
is coupled to supply the signal 5 to the NNAE 36. Based on the signal 5, the NNAE 
36 adjusts the weight data W to adapt the NNAE 36 to generate the pseudo-control 
signal so as to compensate for error A associated with the command control signal 8 C . 

4. Boundedness Statement 

25 The following theorem establishes sufficient conditions for boundedness of the 

error signals and neural network weights in the proposed closed-loop adaptive output 
feedback architecture. 

Theorem 2, Subject to assumptions 1-3, the error signals of the system comprised of 
30 the dynamics in Eq. (7), together with the dynamics associated with the realization of 
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the controller in Eg. (14) and the NN adaptation rule in Eg. (40), are uniformly 
ultimately bounded, provided the following conditions hold 

. Q m >2|C cl ||, X w >c 2 /4, (41) 

where Q m is the minimum eigenvalue of Q. 

5. Processor-Based Embodiment of Adaptive Control System using Direct 

Output Feedback 

Although it is possible to implement the elements 14, 16, 18, 20, and 22 of the 
ACS 10 of Fig. 1 as discrete or grouped analog or digital devices, these elements can 
alternatively be implemented in a processor-based ACS system 10. The processor- 
based system 10 includes a processor 66 and memory 69 storing data and a control 
program, to implement the elements 14, 16, 18, 20, and 22. More specifically, the 
control program can be implemented as software objects or modules that perform the 
functions of the elements 14, 16, 18, 20, 22 as previously described. The data can be 
parameters such as the NN connection weights W, V and/or basis function(s) <|> that 
are updated by the processor 66, as well as temporary data and intermediate 
calculations, commanded state signal levels, plant output signal levels, etc. The ACS 
10 of Fig. 5 can further comprise bus 70 to which the operator interface unit 26, the 
command filter unit 28, the actuator 32, the sensor(s) 34, the processor 66, and the 
memory 68, are coupled. 

In operation, the sensor(s) 34 generate plant output signal y and supply this 
signal to the operator interface unit 26 via the bus 70. The operation interface unit 26 
generates a signal readable or discernible by the operator. If the operator 24 is human, 
the operator interface unit 26 can generate a display or the like based on the plant 
output state signal y. If the operator 24 is a processor or machine, the operator 
interface unit 26 can convert the plant output state signal y into a format usable by the 
operator. The operator 24 if human produces one or more signals through control 
actions applied to a command filter unit 28. For example, in the case of an aircraft, the 
control actions may be applied to control instruments of the aircraft. Alternatively, if 
the operator 24 is a machine, the command signal produced by the operator can be 
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applied to the command filter unit 28. The command filter unit 28 generates the 
commanded output signal y c and the rth derivative of the commanded output signal 
yW. The command filter unit 28 supplies the signals y ci y^ to the processor 66 or to 
5 the memory 68 at a location accessible to the processor 66. The sensor(s) 34 can 
supply the plant output signal y directly to the processor 66, or to the memory 68 at a 
location accessible to the processor 66 via the bus 70. The processor 66 performs the 
functions of the elements 14, 16, 18, 20, 22 to generate a command control signal 8 C . 
The processor 66 is coupled to supply the command control signal 5 C to the actuator(s) 

10 32 via the bus 70. The actuator(s) 32 perform control of the plant 12 in a manner that 
can affect the plant state(s). The sensor(s) 34 sense and generate the plant output 
signal y for the next control cycle. Processing performed by the processor 66 in 
executing its control program can be repeated over successive control cycles as long 
as required to control the plant 12. 

15 6. General Method of the Invention 

Fig. 6 is a flowchart of processing performed by the ACS 10 of Figs. 1-5. In 
step SI of Fig. 6 the method begins. In step S2 a command control signal 5 C is 
generated by inverting an approximate model of the plant dynamics, based on a 
pseudo-control signal v and the plant output signal y. In step S3 the command control 

20 signal 5 C is supplied to control the plant. In step S4 the plant output signal y is 
generated the sensors. In step S5 a tracking error signal y is generated by 
differencing corresponding signal components of the commanded output signal y c and 
optional derivative(s) thereof, and the plant output signal y. In step S6 a pseudo- 
control component signal is generated based on the tracking error signal y using 

25 the transfer function N dc (s)/D dc (s) . In step S7 a transformed signal (Note change 
here) y ac j is generated based on the tracking error signal y using transfer function 
N a d(s) /Dd c (s). In step S8 the rth derivative of the commanded output signal y** 1 * is 
generated. In step S9 the signal y ad is generated to render the transfer function from 
the signal v a d to the signal y^ strictly positive real by appropriate choice of N ad (s) 

30 . In step S10 a training signal 8 is generated by filtering basis fimction(s) § and 
multiplying the filtered basis fiinction(s) § by the transformed signal y &d . In step Sll 
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connection weights W of a neural network are updated in a bounded manner using the 
training signal 8. In step S12 delayed versions of the pseudo-control signal v are 
generated. In step S13 delayed versions of the plant output signal y are generated. In 
5 step S14 the adaptive control signal v a d is generated based on the pseudo-control 
signal v, delayed versions of the signal v, plant output signal y y plant output signal 
y d9 connection weights W, V, and basis function(s) <|) updated based on the training 
error signal 5. In step SI 5 a pseudo-control signal v is generated based on the rth 
time-derivative of the commanded output signal y c (r \ pseudo-control component 
10 signal v dc , and adaptive control signal v ad . In step S16 the method of Fig. 6 ends. 

7. Example of Implementation of the Adaptive Control System having Direct 

Output Feedback Control 

15 The performance of the ACS 10 using output feedback is demonstrated by 

considering the following nonlinear system, consisting of a modified Van der Pol 
oscillator coupled to a lightly damped mode 

x,=x 2 (42) 

20 x 2 = -2(x? - l)x 2 - x, + u (43) 

x 3 =x 4 (44) 

x 4 =-x 3 -0.2x 4 +x 1 (45) 

y = x, + x 3 (46) 

25 The output y has a relative degree of r = 2. From a practical perspective, the 

system can be thought of as a second order nonlinear plant model, whose realization 
consists of states Xj and x 2 , in which the output is modeled as y = However, the 
system contains also a very lightly damped unmodeled mode, with a natural frequency 
equal to that of the linearized plant. This mode is excited by the plant displacement 

30 state (x,) and is coupled to the measurement. 

The output signal y does not have a full relative degree in the presence of the 
unmodeled mode. The low natural frequency of this mode is encompassed by the 
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bandwidth of the controlled system. This introduces a challenging control problem, in 
particular for methods that require the output to have a full relative degree. Moreover, 
this example is treated as if even the Van der Pol model is unknown, and only the fact 
5 that r = 2 is given. This is not an unreasonable assumption in that in many systems, 
the number of plant states and hence the value of r can be deduced from knowledge of 
the behavior of the plant. Thus, the controller design is performed assuming y = u , 
implying that in Fig. 1 the plant transfer function from the pseudo-control signal v to 

y is l/s 2 . 

10 A first order lead-lag compensator was selected to stabilize the associated error 

dynamics. In addition, the first design approach described in Section 3.5 was used to 
satisfy the SPR condition. The resulting two outputs of the compensator are given by 



which places the closed loop poles of the error dynamics at - 3,-1 ± j . The low pass 
filter 60 T -1 (s) discussed in Eq. (26) was chosen as 

T - 1 (s) = -^— . ( 
W 0.5s + l 

20 

It is easy to verify that the transfer function G(s)T(s) of Eq. (31) is SPR. 

A Gaussian Radial Basis Function (RBF) NN with only three neurons and a 
bias term was used in the adaptive element. The functional form for each RBF neuron 
25 was defined by 




(50) 
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0 / (7) = ^(^-7c/) r («-«c/)/^ 2 > o- = V2, i = l,2,3. 



(52) 
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The centers r] ci ,i = 1,2,3 were randomly selected over a grid of possible values for the 
vector r\ . All of the NN inputs were normalized using an estimate for their maximum 
values. The current and two delayed values for the plant output signal y and only the 
5 current pseudo-control signal v were used in the input vector to the neurons. The 
complete input vector consisted of these values together with a bias term, as illustrated 
in Fig. 4. Thus, there are a total of four NN weights in Eq. (29). The network gains 
were F= 50 and X w — 1 . 

In the simulation, the initial states of the system were set to Xj(0) = 0.5, 
10 x 2 (0)=2.5, x 3 (0) = x 4 (0) = 0. The system was commanded to follow the output of a 
second order reference model for the MIU 18, designed with a natural frequency of 

i3 CD n =l [rad/sec] and damping ^ = l/ 42 , and driven by a square wave input command 

m signal j/ c . 

lil First, the controlled system performance is evaluated without the unmodeled 

IP 15 mode dynamics, i.e., removing Eqs. (44) and (45) and setting the output y = x P 

,q However, it will be recalled that the controller has been designed given only the fact 

■ « that r = 2. Fig. 7 compares the system response without NN augmentation (solid line) 

| [1 with the reference model output (dashed line), clearly demonstrating the almost 

IH unstable oscillatory behavior caused by the nonlinear elements in the Van del Pol 

j~jf 20 equation. Figs. 8 A and 8B show that with NN augmentation, these oscillations are 

eliminated after a period of about three seconds. This is accounted for by the 
successful identification of the model inversion error by the NN, which is also 
illustrated in Fig. 8B by comparing the NN output (solid line) with the computed 
inversion error (dashed line). 
25 Next, the effect of the unmodelled dynamics is examined. In this case, the 

response without the NN is unstable, and therefore is not shown. The response with 
NN augmentation is presented in Fig 9. It shows only minor performance degradation 
compared to the full relative degree case of Fig. 8A. Since the unmodeled mode is 
well within the bandwidth of the control system (when viewed with v ad = A ), this 
30 demonstrates that the adaptive system learns to interact with the added mode to 
achieve good tracking performance. 
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In Figs. 8A, 8B, 9 the NN based adaptive controller exhibits a steady state 
tracking error. This error can be removed by introducing an additional integral 
control action when designing the linear compensator. The performance of the 
5 controller with integral action is shown in Figs. 10A and 10B. The steady state 
tracking error is zero, while the transient response behavior is only slightly 
compromised. The bounded NN weight time histories are also depicted in this figure, 
showing that most of the weight adaptation takes place when the command reverses 
direction. 

10 8. Additional Considerations 

The stability results are semiglobal in the sense that they are local with respect 
to the domain D. If the NN universally approximates the inversion error over the 
whole space 9t n+1 , then these results become global. 

The NN update laws consist of a modified gradient algorithm along with the 
15 standard a-modification term as described in the Kim and Lewis J 998 publication. 
These laws have been proven to be passive in Lewis, F. (1999), Nonlinear network 
structures for feedback control, Asian Journal of Control, l(4):205-228. 

The NN learning takes place on-line, and no off-line training is required. No 
assumption on persistent excitation is required. 
20 The ultimate bound for the tracking error can be made smaller by increasing 

the linear design gains. This will result in increased interaction with unknown or 
unmodeled plant dynamics. However, Theorem 2 remains valid so long as 
assumptions 2 and 3 hold. 

In the case of plants of unknown dimension but with known relative degree, 
25 the described methodology applies with a slight modification of the input vector to the 
network: a large range of input/output data should be used, i.e., n x » n . 

9. Conclusion 

The described ACS 10 presents an adaptive output feedback control design 
procedure for nonlinear systems, that avoids state estimation. The main advantage is 
30 that the stability analysis permits systems of arbitrary but known relative degree. The 
full dimension of the plant and its internal dynamics may be known or poorly 
modeled. Only mild restrictions regarding observability and smoothness are imposed. 
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Consequently, the result is applicable to adaptive control of nonlinear systems with 
parametric uncertainty and unmodeled dynamics. 

Any trademarks listed herein are the property of their respective owners, and 
5 reference herein to such trademarks is intended only to indicate the source of a 
particular product or service. 

The many features and advantages of the present invention are apparent from 
the detailed specification and it is intended by the appended claim to cover all such 
features and advantages of the described methods and apparatus which follow in the 
10 true scope of the invention. Further, since numerous modifications and changes will 
readily occur to those of ordinary skill in the art, it is not desired to limit the invention 
to the exact implementation and operation illustrated and described. Accordingly, all 
suitable modifications and equivalents may be resorted to as falling within the scope 
of the invention. 
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